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High Temperature Thermodynamics of the Ferromagnetic Kondo-Lattice Model 
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We present a high temperature series expansion for the ferromagnetic Kondo lattice model in the 
large coupling limit, which is used to model CMR perovskites. Our results show the expected cross- 
over to Curie-Weifi behavior at a temperature related to the bandwidth. Estimates for the magnetic 
transition temperatures are in the experimentally observed range. The compressibility shows that 
the high temperature charge excitations can be modeled by spinless fermions. The CMR effect itself, 
however, warrants the inclusion of dynamic effects and cannot be explained by a static calculation. 
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The discovery of colossal magnetoresistance (CMR) in 
doped rare-earth manganites has attracted considerable 
attention |)-||. The double-exchange (DE) model g has 
long been considered appropriate for describing these sys- 
tems. In this model spins are associated with localized 
electrons, which are coupled ferromagnetically to itin- 
erant e g electrons through a large Hunds coupling (Jh)- 
The kinetic energy of these itinerant electrons may be 
lowered by a ferromagnetic alignment of the ti g electron 
spins resulting in an effective ferromagnetic coupling. 

Most studies of the DE model in the large Jh limit be- 
gin by mapping it to a spinless fermion model, where the 
hopping matrix elements depend on the relative orien- 
tation of neighboring spins |^ ^. Then various types of 
mean-field theories are used to solve the resulting many- 
body problem 0-0. It was pointed out by Millis et al. 
H that in order to obtain the experimental value of the 
transition temperature phonons should play an impor- 
tant role. There also exists strong evidence that in the 
insulating phase above T c charge transport is controlled 
by the motion of small magneto-elastic polarons jll} . 

Notwithstanding the importance of a coupling of lat- 
tice, charge and spin degrees of freedom, the magnetic 
properties of the plain DE model are unusual by them- 
selves, which was shown in finite size calculations and 
variational treatments It was found that in 

Id the ground state showed unusual odd-even effects, 
and that the excitation spectrum even for ferromagnetic 
ground states is very unusual JO]. Similar effects were 
found in higher dimensions, which could be understood 
in terms of the degeneracies of the finite-system fermi 
surface |Ti| . 

In view of the surprising results in fll2|-p"4j for the 
ground state and low energy excitations of the DE model, 
it is necessary to reconsider its finite temperature prop- 
erties with as little approximations as possible. We 
have therefore developed a high temperature series ex- 
pansion (HTSE) for the spin-half ferromagnetic Kondo- 
lattice model in the large Jh limit for the internal en- 
ergy, the magnetic susceptibility, the compressibility and 



the magneto-compressiblity. The calculations have been 
done for the simple cubic lattice up to order j3 w and, 
to compare to the simpler mean-field behavior, for the 
Bethe lattice (defined as the interior of a Cayley tree with 
varying coordination number q) up to order /3 14 as a func- 
tion of the electron density < p < 1 . The ferromagnetic 
Kondo-lattice model is defined by the Hamiltonian 
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Here 5, is the local spin coupled via the Hund's rule 
coupling Jh to the itinerant electron spin a a f3c\ a c i „. We 
are restricting ourselves to local spins of length s = 1/2. 

We are interested in studying this model in the Jh — > 
co limit, in which case there are only five basis states 
\s,m > per lattice site. Two of these correspond to the 
up and down spin states of the local spin when the itin- 
erant electron is absent (a spin-i object), and the other 
three correspond to the case where the itinerant electron 
is present and forms a spin-1 object together with the 
local spin. These states are >, § >, 1 1 , 1 >, 
1, > and |1, — 1 > and can be labelled uniquely by their 
m-quantum number alone. In order to consider a doped 
system, we add a term fi'^2 i a rii >cr to (1) and adjust the 
chemical potential ju such that these five states are de- 
generate (analogous to the HTSE of the t— J model 
leaving the hopping matrix element t as the only energy 
scale in the problem. We set t = 1 . The action of the ef- 
fective Hamiltonian in the Jh — * 00 limit can be written 
||, m' >,= -!f[A n 
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where represents the fermion sign arising from in- 
terchanging fermions on sites i and j and A r , 



m')(l - m), B m , m , =yf$ + m')(l + m) |3|. 
HTSEs can be developed for the quantities of inter- 
est by a cluster expansion method |^6|. In the ther- 
modynamic limit, extensive quantities X are written as 
a sum over all topologically distinct graphs g, X = 
N s L(g) x W(g), where L(g) is called the lattice con- 
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FIG. 1. Plot of (3/Ax vs P for P = 0-5 on the Bethe lattice 
with q — 6 and on the s. c. lattice, and the respective PAs. 
The inset shows the second derivative of the respective PAs 
of /3/A\ with respect to /3, i. e. a = 9|(/3/4x)- 



and is given by the relation, W(g) — X g — J2 g 'cg W{g')- 
Here, X g is the appropriate quantity calculated for the 
finite graph and the sum runs over all subgraphs of the 
graph g. For a graph with L bonds, the weight W(g) is 
at least order (3 L . There is an additional symmetry in 
our problem, such that for graphs without closed loops 
the weight of a graph with L bonds is order /3 2L . Thus by 
including all graphs with up to L bonds one can calculate 
the desired quantities for the simple cubic lattice to order 
f3 L and for the Bethe-lattice ( which has no closed loops 
) to order /3 2L . The calculation of the weights is the time 
consuming step in the calculation. The calculations were 
performed on IBM590 workstations and took about 80 
CPU days in total. 

After transforming back to the density as the relevant 
variable we obtain various thermodynamic quantities as 
series in p (the e s -density), and in f3. E. g., the reduced 
susceptibility can be written as 

4X//3 = 1 + \ P + P(l ~ P) £ ( J2 a ltm p l ) /3 2m (2) 

' m \i=0 / 

The coefficients for this series and for the other quantities 
of interest are available on request. 

To locate the expected paramagnetic to ferromagnetic 
transition temperature T c let us first consider the sus- 
ceptibility series. Although the series for x/P contains 
only even powers of /?, we observe the expected recov- 
ery of Curie- WeiB behavior below a typical temperature 



linearity of /3/x vs. (3 shown in Fig. 1 for both the sim- 
ple cubic (s. c.) and the Bethe lattice with q — 6 for 
p = 1/2 (Similar behavior also exists for all < p < 1). 
Unfortunately related to this crossover from quadratic to 
linear behavior an unphysical singularity appears on the 
imaginary axis. For the series shown in in Fig. 1 this sin- 
gularity is at /3| =6 ~ ±.94i, (3 S S C - w ±.85« obtained from 
a ratio analysis of the respective series. This singularity 
dominates the analytic behavior of the series and renders 
conventional methods for extracting T c , like dlog Pades 
or Neville tables, unsuccessful. Pade approximants (PAs) 
actually allow for an extrapolation of the series beyond 
its normal radius of convergence (see Fig. 1). 

For the g-coordinated Bethe lattices we can use the ex- 
pected mean-field critical behavior — the susceptbility 
of the ferromagnetic Heisenberg model for all T > T c is 
X = const./ (T — T c ) — to obtain T c from a fit of the PAs 
of x^ 1 to a linear behavior in T (see Fig. 1). There is, of 
course, some ambiguity in identifying the region in tem- 
perature where linear behavior holds. We use the lowest 
possible range in T, where the two highest order PAs 
agree to obtain an estimate for T c as a function of p and 
coordination number q. For p = 1/2, this fit gives e. g. 
T« =6 = .196 ± .02, T?= 32 = .485 ± .023, T« =256 = 1.35 ± 
.04 (The errorbars are subjective estimates), which in- 
deed follows the expected scaling with q, i. e. T c 9 oc q 1 ^ 2 , 
reasonably well, and in retrospect validates the PAs and 
our analysis. Since for the Bethe lattices we have enough 
terms, we can actually also use integral PAs |L7|]. This 
yields e. g. Tj? =6 — 0.21 for p — 1/2, which agrees within 
the errorbars with the linear fit. It is worth pointing out 
that even the highest estimate for T c from the integral 
PAs, T c — 0.22i, is very low compared to values obtained 
by Millis et al. ||. If one takes values of the bandwidth 
between 1-2. 5eU one obtains estimates for T c between 
180JT and 460-R", which are well within the experimen- 
tally observed range. 

The analysis of the series on the simple cubic lattice is 
more difficult, since one expects a non-trivial suscepti- 
bility exponent 7 > 1, and also because the series is too 
short to allow for a sensible analysis using integral PAs. 
However, just by comparing the susceptibilities for the 
q — 6-Bethe lattice and the s.c. (see Fig. 1) one notices 
that an eventual intercept with the axis would happen at 
a larger f3 value for the s. c. lattice, resulting in an even 
lower T c . 

If one assumes that the critical region is small then 
one can analyze the series for the s.c. lattice in a similar 
way as before. However, the resulting estimates for T c 
are only rough estimates and provide lower bounds (due 
to 7 > 1). The values at which the PAs still agree can 
certainly be taken as upper bounds. Estimates for T c (p) 
are shown in Fig. 2. The values are consistently about 
25% lower than the q = 6 Bethe lattice estimates. 

The small curvature in the inset of Fig. 1 may be 
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FIG. 2. Plot of estimates of T c (p) for the s.c. lattice 
obtained from various fit intervals as described in the text. 
The spikes arise from instabilities in the PAs, and have no 
physical significance. The thin lines are scaled plots of p(l — p) 
(dash-dot) and (p(l - p)) 1/2 (solid). 
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FIG. 3. The inverse compressibility (p 2 K) _1 as a function 
of temperature T. Shown are the [4,6] -PAs to the 10 th order 
series on the s. c. lattice. The thin lines through each pair 
depict the values of (p 2 ft) -1 from the HTSE for free spinless 
fermions on a s.c. lattice at the corresponding densities p. 



to assume that at least at i -filling the system consists 
of ferromagnetically coupled spin-is and spin-Is, and if 
one assumes some kind of charge ordering (e.g. in an AB- 
structure) this would give rise to behavior reminiscent of 
a ferromagnetic "ferrimagnef . For such a system 

4 X = [(1 + b - P )T + 2^-p{\ -p)T c ]/(T 2 - T c 2 ) (3) 

within mean-field theory |jjgj| . The density dependence 
of (||) is obtained in a mean-field picture. We imagine 
that large spins of size 2S\ = (1 — p) (the empty sites) 
and size 2S*2 = 8p/3 (the occupied sites) are coupled 
via a ferromagnetic J12. Possible other couplings, like 
Jn and J22 are zero in our case of the DE limit, i. c. 
there are no induced spin-spin interactions between two 
empty or between two occupied sites. For p = 1/2 the 
slope of x actually agrees with (||) (for the Bethe lat- 
tices). However, the actual difference from the value for 
a Heisenberg model are small, and for other values of 
p the agreement of our series with (|jj|) becomes worse, 
although still showing the trend contained in (J3J) . Ferri- 
magnetic mean-field theory also predicts the concentra- 
tion dependence of T c — Ji 2 y/6p(l — p)/9 as opposed to 
the expected T c (p) oc p(l — p), see e. g. P,^9{. Apart 
from the unphysical spike around p = 0.3 the agreement 
of estimation from our HTSE with the prediction from 
ferrimagnetic mean-field theory is remarkable. 

This still leaves open the question of emerging struc- 



ing at this within our framework is to consider the com- 
pressibility defined as k = (f^)/p 2 and shown in Fig. 
3. Again we find a quasi-linear behavior below a certain 
temperature. By comparing with a HTSE for spinlesss 
fermions on the s. c. lattice we see that differences in the 
compressibility between the DE model and noninteract- 
ing spinless fermions are small and only become relevant 
at low temperatures. This means that the DE model at 
high temperatures behaves pretty much like a system of 
non-interacting spinless fermions at least with respect 
to charge excitations. We do not see any evidence for 
unusual behavior in k that would be indicative of some 
kind of charge ordering. However, this is not conclusive 
since one would expect unusual behavior to occur at a 
specific wavevector related to the filling, and hence it is 
not surprising to see nothing in our q = compressibility. 

To connect more closely to the CMR effect it would 
be necessary to calculate the magnetoresistance (MR) 
p c (h) j p c {h = 0) — 1, where p c is the resistivity and h the 
applied magnetic field. The MR arises from two sources. 
One is a change in the scattering time due to the mag- 
netic field, which modifies the dynamical properties, the 
other is the dependence of the bandwidth on the mag- 
netic field, which is a static effect. The latter is large in 
the DE model as the effective bandwidth is strongly de- 
pendent on the alignment of the local spins. In a fermi- 
gas, the compressibility is proportional to the effective 
mass. Thus, we can capture this bandwidth effect in the 
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FIG. 4. The MC as defined in the text as a function of j3. 
Shown are two PAs (upper: [4,6] lower: [6,4]) for the respective 
densities p. 



timates for the magnetic transition temperature are close 
to the experimentally observed ones. Mean-field theory, 
as mimicked by the Bethe lattice series, overestimates 
T c by about 25%. There are also some unusual features 
associated with the charge degrees of freedom as indi- 
cated by the slightly ferrimagnetic behavior, but our cal- 
culations are not sufficient to study long-range charge 
order. Nevertheless, the tendencies towards charge or- 
dering should imply a large sensitivity of the CMR ma- 
terials towards Jahn-Tcllcr ordering effects. In further 
work we will investigate the wave-vector dependence of 
the compressibility to address this question. Our results 
for the magneto-compressibility indicate that dynamic 
effects are responsible for the CMR effect, and that a 
static description may be inadequate. 
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magneto-compressibility (MC) K(h)/n(h = 0) — 1. In 
Fig. 4, we show the MC for various densities as a func- 
tion of inverse temperature. As expected, and desired, it 
is always negative. For intermediate p we observe satu- 
ration of the MC around (3 = 1 corresponding to the ap- 
pearance of linear behavior of This does not agree 
with the behavior observed in the MR of the real mate- 
rials. Also the concentration dependence of the MC is 
inconsistent with that of the MR; the higher the concen- 
tration of e g electrons the larger is the MC. This can be 
understood by noting that for larger p there is effectively 
more spin to be polarized by the magnetic field. One way 
to resolve this discrepancy would be to assume that the 
dynamic effects in the MR are more important than the 
static ones captured within the MC. 

In this letter we have provided the first exact results 
for the ferromagnetic Kondo lattice model in terms of a 
HTSE. Although we have used a smaller spin value than 
is realized in the CMR materials, the ti g spin length is 
set to 1/2 instead of 3/2, this should not affect T c much 
as, to leading order for Jh — > oo, the spin-length drops 
out of the problem. The analytic structure of the re- 
sulting series is governed by a singularity on the imag- 
inary axis, which we indentify as being responsible for 
the crossover from the high-T range to the physically in- 
teresting temperature range below the bandwidth. This 
behavior is probably also typical of other itinerant sys- 
tems, although the DE model in the Jh — > oo limit is 
unusual, because the local moment JU < (^i) 2 > is a 
constant, and does not depend on temperature. Our es- 
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